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Using a three-dimensional formalism that includes relativistic kinematics, the effects of negative- 
energy states, approximate boosts of the two-body system, and current conservation, we calculate 
the electromagnetic form factors of the deuteron up to Q 2 of 4 GeV 2 . This is done using a dynamical 
boost for two-body systems with spin. We first compute form factors in impulse approximation, 
but then also add an isoscalar meson-exchange current of pion range that involves the "/n contact 
operator associated with pseudovector -nN coupling. We also consider effects of the p-n'y meson- 
exchange current. The experimentally measured quantities A, B, and £20 are calculated over the 
kinematic range probed in recent Jefferson Laboratory experiments. The pn-y meson-exchange 
current provides significant strength in A at large Q 2 and the jtt contact-term exchange current 
shifts t2o, providing good agreement with the JLab data. Relativistic effects and the yir meson 
exchange current do not provide an explanation of the B observable, but the pity current could 
help to provide agreement if a nonstandard value is used for the tensor pN coupling that enters this 
contribution. 



I. INTRODUCTION 

The deuteron is a bound state of two spin-1/2 parti- 
cles. Consequently, the combined dictates of Poincare 
invariance, current conservation, and parity tell us that 
it has three independent form factors. Conventionally 
these are taken to be the deuteron charge, magnetic, 
and quadrupole form factors, Gc, Gm, and Gq, which 
are related to the Breit-frame matrix elements of the 
deuteron electromagnetic current operator A.n in the 



three deuteron magnetic sub-states 
via the formulae: 
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(Here rj = Q 2 /(4M|), and Q 2 = -q 2 is the absolute 
value of the square of the four-momentum transfer to 
the deuteron.) Three experimental quantities are there- 
fore required to disentangle the electromagnetic current 
of this nucleus. Two of these — the structure functions 
A and B — can be obtained from the electron-deuteron 
differential cross-section using the usual Rosenbluth sep- 
aration: 



da da 
dfl dfl Mott 



A(Q 2 ) + B(Q 2 ) tan 2 
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The structure functions A and B are related to Gc, Gq, 
and Gm, as follows: 



A = Gl + ^Gl + ^G 2 ,, 



(5) 



B = 3,(1 



V)G 2 M . 
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The third observable of choice is the dependence of 
the scattering on the (tensor) polarization of deuterium. 
This can be measured either using a polarized deuteron 
target (and unpolarized beam), or with an unpolar- 
ized target by measuring the polarization of the recoil 
deuterons. Both types of experiment result in the same 
tensor-polarization observable, and in particular both 
can measure: 
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Recent experiments at the Thomas Jefferson National 
Accelerator Facility (JLab) have probed the electromag- 
netic form factors of the deuteron at large space-like mo- 
mentum transfers, too has been measured at Q 2 up to 
almost 2 GeV 2 0, |2L B out to about 1.3 GeV 2 , and 
A to Q 2 = 6 GeV 2 0, 3 . Two recent reviews provide 
up-to-date information on these experiments and their 
theoretical interpretation 0,0]. 

The kinematic range of these data pushes the limit 
of theoretical descriptions of this simplest of nuclei. At 
small values of Q 2 the appropriate degrees of freedom for 
this description are nonrelativistic nucleons, interacting 
via static nucleon-nucleon interactions, with small cor- 
rections due to meson-exchange currents and relativistic 
effects 0, H H E EH Provided Q 2 is below the scale 
of chiral symmetry breaking, contributions to electron- 
deuteron scattering can be organized using the nuclear 



2 



effective theory first proposed by Weinberg 0, E, 0] , 
and the results are in good agreement with the extant 
data below Q 2 ~ 0.5 GeV 2 . As Q 2 increases to 
~ 1 GeV 2 , nucleon and meson degrees of freedom may 
still be appropriate but a relativistic formalism is needed 
to account for relativistic kinematics and dynamics. At 
very high values of Q 2 one expects a transition to a 
regime where quark and gluon degrees of freedom provide 
the most natural description of data. To date, the ex- 
perimental and theoretical situation for elastic electron- 
deuteron scattering can be summarized by saying that 
there is no significant evidence of having reached the 
quark regime |5j. The most important aspects of the 
internal quark structure of a nucleon appear to be taken 
into account via nucleon form factors and relativistic NN 
dynamics. 

Considerable effort has been invested in the construc- 
tion of these relativistic formalisms for the NN bound 
state and associated meson-exchange currents. The goal 
is to address elastic electron-deuteron scatter ing at Q 2 o f 
a few GeV 2 , see e.g. Refs. [H E3 E E H H 
Such an approach is a logical extension of the stan- 
dard nonrelativistic treatment of the NN system (see, 
e.g. Ref. |24[) which has had a significant amount of suc- 
cess in describing A, B, and £20 using a non-relativistic 
NN interaction that is fit to the NN scattering data 
(see, e.g. Refs. 0, 1 0). Like the non-relativistic ap- 
proach, relativistic approaches to the NN problem are 
grounded in a phenomenological description of the NN 
scattering data. However, they seek to implement dy- 
namics that obeys the Poincare algebra — even if only 
approximately — and in so doing they go beyond the non- 
relativistic treatment of the NN system. If electro- 
magnetic interactions with the deuteron are also to be 
considered then — regardless of the momentum transfer 
involved — it is crucial that the consequences of electro- 
magnetic gauge invariance be incorporated in the calcu- 
lation. Minimally this means that the electromagnetic 
current of the deuteron must be conserved. Indeed, the 
derivation of Eqs. Q-© assumed that the deuteronic 
current A^ was conserved. This motivates the use of 
meson-exchange interactions for which methods of quan- 
tum field theory may be used to construct the conserved 
current. 

The three-dimensional "equal-time" (ET) formalism 
which was developed and applied in Refs. 0, 0, 
is one such method. This approach starts from the four- 
dimensional Bethe-Salpeter formalism and the Mandel- 
stam construction of the electromagnetic current. It in- 
cludes relativistic kinematics, negative-energy states and 
relativistic pieces of the electromagnetic current explic- 
itly at all stages of the calculation. In Ref. 0] we 
reported on impulse-approximation calculations of the 
form factors The Bonn-B interaction [25|, that 

was fit to NN scattering data using a relativistic equa- 
tion with only positive-energy states, was used. In con- 
trast to Ref. [23 we employed pseudovector ttN coupling, 
since it is then easier to implement the constraints of chi- 



ral symmetry. In Ref. [18| we focused particularly on the 
role of negative-energy components of the deuteron wave 
function in the ET formalism. The NN interaction with 
negative-energy components was obtained by adjusting 
the aNN coupling so that the deuteron binding energy 
was the same as in Ref. [25|. Once this was done the 
inclusion of negative-energy components of spinors pro- 
duced only modest effects: there were noticeable changes 
to some deuteron observables at larger momentum trans- 
fers but not much improvement in the description of the 
experimental data. This is i n ag reement with earlier re- 
sults of Hummel and Tjon 19], and is not surprising 
since pseudovector coupling suppresses the coupling to 
the negative-energy states. 

The results of Ref. for observable £20 were in 
good agreement with experiments but the description 
of the deuteron magnetic form factor was poor. Sim- 
ilar conclusions held regardless of whether or not the 
negative-energy components of the deuteron wave func- 
tion were included. Some other relativistic approaches, 
e.g. Ref. |23|, have found larger effects for negative-energy 
components and this point is discussed further in Sec. lIXI 

However, a number of important effects were not 
included in the impulse-approximation calculation of 
Ref. (hereafter "ETIA"). These included the dynam- 
ical boost of bound-state vertex functions, an isoscalar 
pionic meson-exchange contribution to A^ that arises 
from use of pseudovector nN coupling, and the pnj 
meson-exchange current. In the present work we use the 
positive-energy-state ETIA calculation as a baseline. We 
then include the various effects listed above that were 
omitted in Ref. and also display the contribution of 
Z-graphs computed to first order in perturbation theory. 

These different contributions to A^ can be organized 
accordin g: to the power counting developed by Wein- 
berg [T2I A1$L Hif that has already been successfully ap- 
plied to a number of electromagnetic reactions involving 
deuterium. Since the power counting is based on an ex- 
pansion in powers of 

v = IpI,™.*, \g\ ^ 

with A the scale of chiral symmetry breaking, it breaks 
down at Q 2 ~ 0.5 GeV 2 , but examining the different 
contributions to A^ using the nuclear effective theory 
provides a way to anchor their low-momentum behavior 
in a systematic way. 

In this counting the calculation of Ref. included 
all mechanisms of O(e) and 0(eP). Formally the most 
important neglected mechanism occurs at 0(eP 2 ), and 
arises because electromagnetic currents are evaluated in 
the Breit frame, where the deuteron has momentum — iq 
in the initial state and momentum in the final state. 
In this case the deuteron wave functions must be boosted 
from the deuteron rest frame to obtain the initial and 
final-state wave functions in the Breit frame. The ET 
formalism developed in Ref. [Trl ITtI Il8| was not boost in- 
variant, and this failure to respect Poincare invariance re- 
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suited in an error at 0(eP 2 ) in the calculation of Ref. . 
In order to remedy this, a boost rule for scalar particles 
was developed in Ref. (2fi| where it was shown to provide 
an exact result for the boost of the two-body energy. 
In this work, the boost rule is extended approximately 
to the case of two interacting spin- 1/2 particles as dis- 
cussed in Sec. We have verified numerically that the 
ET equation together with the boost rule for the NN 
interaction provides eigenstates of the deuteron with en- 
ergy Ep = V M 2 + P 2 when the total momentum is P. 
The deuteron mass M is therefore invariant for the total 
momentum values of interest in this work. 

An exact solution of the boost problem has been given 
in the instant form of dynamics by Coester and Poly- 
zou 1^3, following the methods suggested by Bakamjian 
and Thomas |28[ . In that solution, interactions are in- 
troduced solely into the rest-frame mass operator. In 
contrast, here we use the ET reduction of field-theoretic 
expressions to find an effective hamiltonian as well as 
currents that are consistent with that hamiltonian. (In 
particular, we ensure that the currents satisfy a Ward- 
Takahashi identity.) In the two-body rest-frame system, 
the hamiltonian found via our ET reduction could be 
used to define a suitable mass operator, which could then 
be used to implement an exact boost to other frames fol- 
lowing the methods of Ref. [23 ■ We do not employ this 
"exact boost" construction because the electromagnetic 
currents that are obtained from the ET reduction are not 
consistent with it. Instead, in our work we develop an ap- 
proximate boost rule which maintains a clear correspon- 
dence to field theory, and so facilitates a more consistent 
treatment of strong and electromagnetic interactions. 

Another important contribution not included in 
Ref. enters at 0(eV 3 ). When pseudovector (PV) 
7rV coupling is employed an additional contribution to 
the two-body isoscalar charge operator — first identified 
by Riska [2^ — must be included. Without it the unitary 
equivalence between PV and pseudoscalar irN coupling 
will not be respected |30| . In this work we include this 
PV-coupling current and find that it is significant, albeit 
not as important as when it is added in some nonrcla- 
tivistic approaches • The decreased importance of this 
effect in our approach is traceable to a factor of half in 
the PV-coupling current in the ET formalism relative to 
the current used in Ref. pj. This factor must be taken 
into account if the unitary equivalence between differ- 
ent techniques for obtaining relativistic corrections to the 
AW interaction and currents is to be maintained |31| . 

The pirj meson-exchange current (MEC) has often 
been invoked in calculations of electron-deuteron scat- 
tering. In this MEC the photon interacts with the p- 
meson cloud of one nucleon, producing a pion which is 
absorbed by the other nucleon. Thus the 7A interaction 
is of short range but the MEC has long range because 
a pion exchange is involved. This MEC contributes to 
Gm at 0(eP 4 ) in the nuclear effective theory, but as we 
show it can produce substantial effects in the magnetic 
form factor at large Q. We include it in our relativis- 



tic calculation and consider the extent to which it can 
help to explain the magnetic form factor in the ET for- 
malism. Although there is not much improvement in the 
description of data if the standard tensor pN coupling of 
one-boson-exchange JVAf models is adopted, other values 
of f p /g p may help to explain B(Q 2 ). Note that we do 
not include the 0^177 current that has been considered by 
Hummel and Tjon [19( because it is of shorter range. 

Our paper is structured as follows. First we give a brief 
review of the ET formalism in which we display expres- 
sions for the bound-state equation, the NN interaction, 
and the current matrix element — all for the particular 
case of an instantaneous two-body interaction. Second, 
we discuss the boost of the bound state wave function 
that is needed to evaluate matrix elements in the Breit 
frame. Third, we discuss inclusion of meson-exchange- 
current contributions, especially the PV-coupling current 
and the pirj MEC, each of which has been found to give 
significant contributions to electron-deuteron scattering. 
We then also describe the inclusion of Z-graph effects 
(coupling to the negative-energy states) in perturbation 
theory. Lastly, we present our results for A, B, and £20 
when each of these effects is added to the baseline ETIA 
calculation. 



II. THE EQUAL-TIME APPROACH 

A number of alternative 3D relativistic treatments of 
deuteron dynamics exist (see, for instance, Refs. [l9l I20L 
I2TL I22I l23|h Of these, the formalism that is closest to 
this work is that of Hummel and Tjon [l^, although here 
we eliminate some approximations made in Ref. |l9j |. In 
this section we summarize the development of the ET 
formalism with respect to obtaining the one-body limit 
and a systematic reduction from 4D to 3D. 

Consider the 4D Bethe-Salpeter equation (BSE) for a 
bound-state vertex function, T: 

T = KGqT. (10) 

Here K is, in principle, the sum of all two-particle- 
irreducible NN — » NN graphs. The NN propagator 
Go is the product of spin-half Fcynman propagators for 
each nucleon: Go = id\d,2- In studies of this equation for 
the deuteron bound state the kernel K included a 
set of single-boson exchanges — in analogy to many non- 
relativistic potential models — yielding the "ladder" ap- 
proximation. However, it is well known that in such an 
approximation the Bethe-Salpeter equation does not give 
the correct one-body limit [33j. In other words, if we con- 
sider unequal-mass particles, and take one of them to be 
very heavy, Eq. (|10l) does not become the Dirac equa- 
tion for the light particle moving in the static field of the 
heavy one. This limit is only properly treated in Eq. 1|10|) 
if the full set of ladder and crossed-ladder graphs is taken 
for K 33]. In Ref. 0] we provided a remedy to this flaw, 
and showed that the pieces of the graphs that appear in 
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K and are responsible for the one-body limit can be re- 
summed so that Eq. IjlUI) becomes: 



T = U(Go + Gc)T, 



(11) 



where the precise form of Gc was derived in 0, . For 
exact correspondence between ()l(Jfl and (lllf) we should 
have: 



K = U + UG C K. 



(12) 



At the level of the one-boson-exchange interaction, where 
K and U have only their lowest-order pieces, we see that 
Eq. defines an improved "ladder" Bethe-Salpeter 

equation, which does have the correct one-body limit: 



T = AT (2) GI\ 



(13) 



where G = Gq + Gc- Note that if the field theory to 
be solved involves nucleons and mesons then the approx- 
imation K — > K 1 - 2 ^ is equivalent to restricting ourselves 
to a one-boson-exchange kernel. 

This equation is still four-dimensional. The next step 
is to perform a systematic reduction to three dimen- 
sions. We can motivate the reduction scheme by following 
Salpeter and assuming that the dominant interaction 
is the instantaneous part, i.e., the part obtained by the 
replacement: 



rir 2 



rir 2 



(14) 

where 1^ and T 2 are Dirac-matrix structures of the 
meson- nucleon vertices and q = (qa,q) is the four- 

(2) 

momentum of the meson. Since A~ inst depends only on 
the three-vector q this replacement reduces Eq. I|13l) to 
a three-dimensional equation. If we were to stop at this 
point, the reduction would not be systematic nor would 
it include retardation effects. A central result of Ref. 
was that the three-dimensional reduction could be im- 
plemented in a systematically improvable way. Instead 
of simply neglecting the retardation effects, they can be 
reorganized into a 3D interaction kernel K^t- This in- 

(2) 

teraction shares the property of K^L that it does not 
depend upon the time-component of momentum trans- 
fer. If the reduction includes only positive-energy matrix 
elements, then K^t is the 3D interaction obtained from 
time-ordered perturbation theory. 

Thus the reduction to three dimensions produces the 
following two-body equation: 



Tet — ^et(G)Fet- 



(15) 



Here the three-dimensional propagator (G) is obtained 
by integrating over the time-component of relative four- 
momentum: 



<G)^ J ^[G (p;P) + Gc(p;P)} 



(16) 



and the ET interaction is defined in lowest order by 

(G)K ET (G) = (GKG) . (17) 

(Hereafter we always denote integration over zeroth com- 
ponents of relative four- momenta by angled brackets.) 

In order to explain the significance of the Gc term in 
(G) and its connection to Z graphs and the one-body 
limit, consider the more standard equal-time Green's 
function j^, Eij that omits Gc ■ It is sufficient to work in 
the cm. frame of two particles of equal masses because 
the boost discussed in Sec. will provide the results 
needed in other frames. Then one finds 



<Go> = 



A+A+ A7A 2 



E-2e E + 2e' 



(18) 



where A^ are related to projection operators onto posi- 
tive and negative-energy states of the Dirac equation, E 
is the total energy, and e = (p 2 + m 2 ) 1 / 2 . The propaga- 
tor (Go) is not invertible |37|. since it has no components 
in the H — and — h sectors. This is related to the lack 
of a correct one-body limit in the ladder BSE. If we had 
applied the 3D reduction l|14fl to Eq. (|10fl we would have 
obtained the Salpeter equation: 



Ts = -^inst(Go)rg 



(19) 



which has the non-invertible (Go) in the intermediate 
state. However, adding Gc — which comes from resum- 
ming pieces of the crossed-ladder graphs — before reducing 
to three dimensions gives a 3D JVJV propagator: 



(G) 



Aj-A+ 
E-2e 



2e 



ATAJ 
2e 



A-i A 2 
E + 2e' 



(20) 



This is the three-dimensional propagator that was de- 
rived by Mandelzweig and Wallace [33, here specialized 
to the cm. frame for two equal-mass particles. In its 
more general form for unequal masses, the contribution 
(Gc) to (G) provides the correct one-body limit as either 
particle's mass tends to infinity. The propagator also is 
invertible. 

With regard to Z-graphs, one may compare the +H > 

++ piece of 



K^(G)Kl 2) 



ET 



(21) 



with the amplitude obtained at fourth order in the full 4D 
field theory. We find that the contribution of negative- 
energy states agrees at leading order in 1/m [l8|, and 
that this would be true even were different mass particles 
considered and the mass of either one taken to infinity. 
In other words, effects such as Fig. ^ are included in a 
bound-state calculation that employs Eq. (|15|) . And this 

(2) 

is true even if only the instantaneous ladder kernel L 
is used, because of our careful treatment of the one-body 
limit. 
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shell to determine the value of p\ — p'x°- The interaction 
(|24[l agrees with the p = 1 interaction used by Adam, 
Goller, and Arenhovel [3l|], up to terms of order p 4 /m 4 . 



FIG. 1: One example of a Z-graph which is included in our 
3D equation 1151 . 



III. THE INSTANTANEOUS AW 
INTERACTION 

In Ref. |Tsj we examined the importance of retardation 
effects in K^t- We found that they had little impact 
on deuteron electromagnetic form factors. Hence in this 
work we report only on results obtained using Eq. 115|) 
with an instantaneous interaction. This is consistent with 
use of the Bonn-B interaction that takes an instantaneous 
form in the center-of-mass frame. Other reductions to 
three dimensions are certainly possible but the result that 
retardation effects are small should not be sensitive to the 
reduction used. 

As an example of how the reduction to an instan- 
taneous, three-dimensional interaction is performed we 
consider the four-dimensional pscudovector pion kernel: 



gin ■ t 2 75*71 • (pi - K) 7I72 ■ O2 - p' 2 ) 



9 9 

q z - m* 



2m 



2m 



(22) 



We use the identity: 



u lTx 



5 7i ' (Pi "Pi) 



2m 



(P? 



Pi 



ei 



2m 



•4) 5 . 
7i7i 



?}«!, (23) 



where u'x and u\ are Dirac spinors for positive-energy 
states, and carry out integrations over time-components 
of momenta as implied by the brackets in Eq. I|17l) . To de- 
rive the instantaneous interaction we then take the static 
limit, i.e. assume \E—e—e'\ <C u>. In the language of effec- 
tive field theory this means we consider only the effects of 
"potential pions" (which have go ~ 0) and ignore the im- 
pact of "radiation pions" (qo ~ m^) on our results |39j| . 
The result for the instantaneous ET interaction in 
states is then: 



1ST! 

i\ ET 



9l 



q z + m 
dx - d[ 



7! 5 7l 



-j-^-7i7i72 72 
4m z 



4m 



7i 5 7?72 5 +7i 5 7 2 5 72°)}"i«2, (24) 



where u\ and u 2 are Dirac spinors, dx — Ep — £i(pi) — 
e 2 (p2) and d[ = Ep — ei(p' 1 ) — e 2 (p' 2 ). This interaction 
differs from a pseudoscalar one because of the terms in- 
volving di and d[, which are sub-leading in 1/m. They 
arise because — although we ignored retardation in the 
denominator of Eq. Il22[l — we must still choose an energy 



IV. CURRENT CONSERVATION IN THE 
EQUAL-TIME APPROACH 

In order to compute deuteron electromagnetic form 
factors we must also construct a conserved electromag- 
netic deuteron current. In four dimensions there are two 
pieces to the impulse current, the first being determined 
by the Go part of the propagator: 

f(P + q)G' J ^G T(P) = 

it{P + q)d 1 { Pl )d 2 {p 2 +q)j^d 2 (p 2 )T{P) + (1 <-> 2), 

(25) 

where T is the solution of Eq. (|10|) for initial deuteron mo- 
mentum P = pi +P2, r" is the solution for final deuteron 
momentum P + q, d n is the Dirac propagator for particle 
n, and jp™ is the usual one-nucleon current for particle 
m 



An) 



» 



F 



W(Q2) 7 w+if>)(Q^^y 



(26) 



where = v(l + r 3 ) i s the charge. Using the Ward- 
Takahashi identity associated with formally modified but 
practically identical form of this current |4jj it is easy to 
show that the current I125H is conserved, i.e. 

g"f(P + q)G' Q J ^G T(P) = 0, (27) 

provided that T obeys the ladder Bethe-Salpeter equa- 
tion. 

To obtain an equivalent three-dimensional current, we 
can replace T — > T$, yielding the conserved current 

As* = f' s {(G'^GoK^Go) + (G'oKW'G'oJ^Go) 

-(G J M G )}r s (28) 

where V , G' and K^' include the momentum Q from 
the photon absorption. Since we are considering the 
equal-time case we can also replace K 1 - 2 ^ — ► Ki nst , and 
because i^i ns t has no dependence on time components of 
momenta, that replacement applied to Eq. (|28|l leads to 

•As,n = F' s {(G' J fl G }K iast (Go) 



(G'q) K- ms t (Gq J p, Go ) 

(G' J^Go)}T s (29) 



where the angle brackets indicate where integrations have 
been carried out over time components of momenta. This 
latter expression collapses to the impulse approximation 
form when Eq. (|19|l is used in the first two terms: 



«4s,/j = T' s (G' Q J^Go)Ts- 



(30) 
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The three terms in Eq. I|28(l are equal and they simplify 
to just one term in Eq. (|3U1) . 

However, in this work we did not begin with the ladder 
BSE. Instead we began with the 4D equation (|llfl . Con- 
structing a conserved impulse approximation current for 
the vertex function which is the solution of Eq. (|11|) is a 
little more involved. In Ref. we showed how to add 
a piece to the current l|25|) which results in a conserved 
current when T is the solution of Eq. giving a total 
4D current 

f (P + q)G' J 0<fl G T(P) + f(P + q)G c ^T(P). (31) 

The explicit expression for Gc„ can be found in Ref. [3. 

With this four-dimensional current in hand, we cal- 
culate the reduction to three dimensions as in Eq. (pTTTf> . 
Replacing T by Tet we end up with 

Amst,n = T ET Qi nst ^ T E t (32) 

where: 

£Lt,„ = (G'W = (G' J,G Q + G Cfi ). (33) 

This is analogous to the reduction employed for the 
bound-state equation itself. Once again, this reduction 
can be performed in a systematic fashion, but here we 
keep only the results for an instantaneous interaction. In 
that case Ai ns t,n is conserved, provided that Tet is the 
solution of Eq. 1)1 5H using the instant form of the inter- 
action. The explicit form of G[ nst is 18]: 

(Pi . P2 5 P Q) = * (di (Pi ) d 2 (pa + Q)jW d 2 (pa)) 
+ i(d 1 ( Pl )4(P2 + Q)ji 2 ld c 2 (p 2 )) + (1 <-> 2)(34) 

Only the 7 M piece of is relevant for charge conserva- 
tion, since the piece proportional to a^ v is automatically 
conserved. Meanwhile, d c n is a one-body Dirac propaga- 
tor used in Gc (P) to construct the approximation to the 
crossed-ladder graphs. Correspondingly, d c n appears in 
Gc(P + ?)• It does not equal d£, even if particle n is not 
the nucleon struck by the photon. Finally, 

i® = (»7£> if^^^f, (35) 

e 2 + £ 2 

with p2 = (e(pa)) Pa)' The current defined by Eqs. 
(|35|l includes the effects of photons creating Z-graphs 
by means of couplings from positive-energy states to 
negative-energy states. 

This defines our impulse-approximation current. De- 
tailed results for this current employed with the solutions 
of Eq. I|15fl were presented in Ref. • The various re- 
finements such as the ji 2 ^ term in Eq . (|34[) that were not 
considered by Hummel and Tjon 19] produce only very 
small changes for electron-deuteron scattering. They are 
incorporated in our calculations because they ensure cur- 
rent conservation when the ET propagator is used. 



V. DYNAMICAL BOOST 

To compute the current in Eq. I|32|) we consider elec- 
tromagnetic matrix elements in the Breit frame where 
P = — iq in the initial state and P' = iq in the final 
state. In the impulse approximation, using the instant 
interaction and the instant current in the Breit frame, 
the current matrix element may be expressed as 

f d 3 p - 1 1 „ 1 
J J2nf ^ P + 4 q ' 2 q ^ inst -^ P ~ 2 q ' ~ P; P ' Q > 

r(p-iq;-iq). (36) 

In the ETIA calculation of Ref. [3 Eq. ® is calculated 
using vertex functions Y that are obtained by solving the 
bound-state equation 

{E F -e 1 -e 2 )\E F )=v\E P ), (37) 

where e„ = yjra 1 + p^ , for the state vector \Ep), and 
then computing: 

r(k;P) = (k;P|G 1 | J B P ). (38) 

The operator 

Het = ei + e 2 + v (39) 

may be interpreted as the effective hamiltonian. With 
the usual conditions on the interaction v, it is bounded 
from below and may be used to generate a Hilbert space 
of states. 

The interaction v could be defined only in positive- 
energy states, or it could be the effective interaction in 
positive-energy states that includes the effects of cou- 
plings to negative-energy states: H — , — h and . The 

couplings to negative-energy states obtained from the ET 
reduction do not produce any singularities in the effec- 
tive interaction v. Once the interaction is defined in the 
cm. frame of the two particles, where P = 0, the low- 
est eigenvalue is the rest mass of the two-particle bound 
state, M. 

The problem here is that this does not necessarily guar- 
antee that 

E P = VM 2 + P 2 (40) 

in other frames. That is because solving Eq. I|37|l in an 
arbitrary frame means using the ET interaction K^t for 
v, and Ket is calculated using Dirac spinors with ar- 
guments p ± iP. Such an interaction depends on total 
momentum P, and is not guaranteed to produce a bound 
state of the appropriate energy. In fact, in Ref. 0] we 
found that Ket must be renormalized by a factor A(P 2 ) 
if l|4U|) is to hold in frames other than P = 0. The ad hoc 
factor A varies linearly from 1 at P 2 = to about 1.10 
at P 2 =4 GeV 2 0. 

This occurs because the two-body ET interaction ob- 
tained by solving Eq. I|37|l using the ET interaction does 
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not respect the Poincare algebra. The basic requirement 
of Poincare invariance is that states must transform un- 
der a unitary representation of the Poincare group. The 
ten generators of translations in time, translations in 
space, boosts, and rotations are the hamiltonian oper- 
ator, taken here to be of the form 1|39|) . the total momen- 
tum P, the boost operator K, and the angular momen- 
tum operator J. They obey the following commutation 
relations. 



[P,fl] =0, 
[K,H] = iP, 

[Ji,Pj] = ie ijk P k , 



[J,H}=0, 
[K i ,P j ]=i8 ij H, 
[Ji, Kj] — ie t j k K k , 
[K i ,K j ] = -ie ijk J k . (41) 



The two-body boost problem is to constrain the inter- 
action, v, such that the Poincare generators satisfy the 
commutation rules that are required for Poincare invari- 
ance. Our strategy will be to take v to agree with K^t 
in the two-body rest-frame (P = 0) and then impose 
the conditions (|41|) in order to obtain approximate ex- 
pressions for the matrix elements of v in a frame where 
P =^ 0. We denote the hamiltonian of this "approximate 
boost" procedure by -Bet- 

As mentioned earlier, we could instead employ an ex- 
act solution to the two-body boost problem by following 
Ref. [27j and introducing the interaction K^{P = 0) 
into the rest-frame mass operator M. A hamiltonian 
that satisfies the Poincare algebra would then be: 



H 



BT 



M 2 



(42) 



The reason we choose not to apply Eq. (|42|) is that for 
electromagnetic interactions we also need current oper- 
ators that are consistent with both the rest-frame and 
boosted hamiltonian. The ET formalism solves that 
problem by performing a consistent reduction of the 
quantum field theory. As discussed above, this yields 
currents which are consistent with -Ket- Such a construc- 
tion is significantly more difficult for Hbt- To indicate 
this we merely point out that if -J/et exactly satisfied 
the Poincare algebra then the two different hamiltonians 
i?ET and Hbt would be related by a unitary transforma- 
tion, i.e., there would exist an operator U with U^U = 1 
and: 



i?ET = UHbtU^ . 



(43) 



If this U were known, we could derive currents for use 
with states obtained from the hamiltonian (|42|l by trans- 
forming the currents derived in the ET formalism using 
the transformation l|43[) . However, this unitary transfor- 
mation is not known. So, in order to respect the consis- 
tency between the ET bound states and current opera- 
tors, we will focus on the problem of finding an approx- 
imate boost that preserves the form of Het in frames 
with P ^ 0. 

It should be mentioned that if P 7^ the hamiltonian 
i?ET has small differences from the operator e\ + e 2 + 



-Ket that was used in our analysis of the Ward-Takahashi 
identity. We leave the issue of boosting the currents so 
as to exactly maintain the Ward-Takahashi identity at 
P 7^ as an unsolved problem. 

An approximate solution to the boost problem within 
the subspace of an eigenvalue of the hamiltonian has been 
obtained in Ref. [2(| for the case of two scalar particles. 
It provides an exact result for the energy of the two-body 
system. Here we review the basic elements of that proof 
and extend the result approximately to the case of two 
spin-half particles. 

For the ET (and a number of other 3D) formalism(s), 
the two-body boost problem is the same as for the in- 
stant form of hamiltonian dynamics 41]. In this form 
of dynamics, total momentum and angular momentum 
operators do not depend on the interactions. They are: 



P = Pl +p 2 , 



(44) 



J = rixpx + r 2 x p 2 + -(cti +cr 2 )- (45) 

It follows from H commuting with P and J that v must 
be translationally and rotationally invariant. The com- 
mutation relation of K and H requires the boost opera- 
tor to depend on the interaction. Such a boost is called 
"dynamical" in order to distinguish it from a "kinemati- 
cal" boost for which the boost generator does not involve 
the interaction. The commutation rule between the dy- 
namical boost generator and hamiltonian thus involves 
v 2 terms. 

Bakamjian and Thomas [2^] derived the form of the 
boost operator for instant dynamics. For a free boost, 
the operator is 

K = \ (r l£l + exn) - + (1 - 2). (46) 

2 2 ei + m 

When interactions are present, there is an interaction 
part of the hamiltonian, v, and an interaction part of the 
boost operator that is given approximately by 



K v = - (Rv + vTL) 



(47) 



where R = i(ri + r 2 ). This boost operator omits 
some terms from the standard Newton- Wigner construc- 
tion |22j for R, and likely could be improved. We leave 
such possible improvement of the boost for future analy- 
sis. Although it involves a significant approximation, in 
the case of the deuteron the boost operator of Eq. l|4T|) 
yields the relationship (|4*U|l between the mass and energy 
of the bound state to good accuracy. 

With the definitions given above, [Ki,Pj] = iHSij is 
satisfied. However, there is an error term of order 1 /m 2 in 
the commutation relation [Ki,Kj] = —ieij k J k . Accept- 
ing this error, the interaction v must take an appropriate, 
but unknown, form consistent with [K, H] — iP. 

Noting that the free-boost operator, Ko, and the free 
hamiltonian, Hq = e± + e 2 , obey the commutation re- 
lation [Ko,Hq] = iP, it follows that when interaction- 
dependent terms are introduced into that commutation 
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relation, their contributions must sum to zero, i.e. 



[K ,v} + [K v ,H ] + [K v ,v} = 0. 



This is equivalent to 



where we have defined: 



1 a x x pi 



2 ex + m 
A = ex - e 2 . 

Algebraic manipulations lead to 



+ (1-2); 



(48) 



-[R, H v + vH + v 2 } + -[rA + Ar, v] + [K a , v] = 0, 



(49) 
(50) 



[R, H ]v + v[R, H Q ] + ^H[R, v] + ±[R, v]H 

+ ^(Aru — vrA) + [K CT , v] = 0. (51) 



elements involve 
d 3 p' d 3 p 



(£p|p';P)£ r 



(2tt) 3 (2tt) 3 
{ [A(p'; P) + A(p; P) + B op ] (p' ; P|«|p; P) 

-i[Kp, CT , (p'; P|«|p; P)] - iK(p'- P\v\p; P) 
-H(p';PMp;P)k CT }<p;P|£ P ) = 0,(54) 



where 



A(p;P) 



2E 



i ^ 1 + d^p)> ) w^P) i (55) 



ap 



with d = ei(p;P) + €2(p;P) being the value of H , 
and A(p';P) the same as A(p;P), except with p — > p'. 
Meanwhile, 

d p' P d p P d 

Bop = ap + ^T^ + ^v (56) 



In quantum field theory the boost of a mass eigen- 
state appears to be kinematical. In fact the boost ve- 
locity j3 = P / Ep depends upon the eigenvalue in the 
cm. frame, M, which of course involves the interaction. 
This simple observation provides the key to solving the 
boost problem in instant dynamics. The v 2 terms in the 
dynamical boost may be eliminated when the boost is 
considered within the subspace of a single eigenvalue of 
the mass. 

The boost velocity appropriate to a mass eigenstate 
does not enter the commutation relations. In order to 
restrict the boost to the subspace corresponding to a par- 
ticular eigenvalue of H we evaluate the matrix element of 
(|5T)l between eigenstates of mass M that obey Eq. (|37Jl . 
This yields: 

(E P \ ([R, H ]v + v[K, H ] + ^E P [R, v] + i[R, v]E P 

+ ^(Arv - vrA) + [K a ,v]) \E P ) = 0(52) 



The v 2 pieces in the third and fourth terms of Eq. I|51() 
have been eliminated in favor of the energy eigenvalue 
Ep. Because of this, the boost appears to be kinematical 
in much the same way that the boost of a mass eigenstate 
in quantum field theory appears to be kinematical. 

Performing the manipulations discussed in Ref. |2fi| we 
can convert Eq. (|52|l into an equation that is linear in v, 



{Ep\{\ 



2 v 1 + £) [R> Ho]v + l v[B " Ho] ( 1+ it 

1 1 l r „ , P P P P 

-Ep [R, v] + -[R, v]E P + ^-rv - vr^- 



-^ a ,v]}\E P ) = 0.(53) 



This equation is solved in momentum space in order 
to determine the form of v. Momentum-space matrix 



1 (ax +a 2 ) x P 

2 Ep{E P +M) ' 



(57) 



and 



with 



qi=p { 2-Ep-T¥i U 



(59) 



Derivatives in B op act only on the interaction. 

With spin effects omitted, a solution of Eq. I|54|l was 
obtained such that 



( P ' ; PMp; P) = /(p; P)«(p', p; P)/(p'; P). (60) 



Here: 



(A(p;P)+B op )/(p;P) = 0. 



(61) 



The form of v is deduced from the condition BopS = 
and the boundary condition that, for P = 0, the interac- 
tion must be the cm. frame one. Now constructing the 
rotational scalars 



E 2 

£yp 



Pc'Pc = P P 

one can check that 



(P ; ■ P) 2 

Ep ' 

(P-P)(P / -P) 

Ep ' 



(62) 



B opP 2 = 0, 
B op p1 2 = 0, 

BopPc ■ Pc = 0. 



(63) 
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So, it follows that if v — v c (p' c ,p c ) is an arbitrary func- 
tion of p 2 , p 2 and p c • p c , then 



B op w c (Pc,p c ) = 0, 



(64) 



and the condition B op {j = is satisfied. Thus, in the 



absence of spin, 



where 



and 



w(p',p;P) = v c (p' c ,p c ). 



(p P)P 



(65) 



E P (Ep +M)' 



(p' P)P 

E F (E P + M) ' 



(67) 



In the cm. frame, p c and p' c are the standard relative 
momenta. When the total momentum is in the z direc- 
tion, the ^-component of relative momentum p c is con- 
tracted according to p cz = p z /j, where 7 = Ep/M. The 
components of relative momenta perpendicular to the to- 
tal momentum are unaffected: p c ^ = p^ . The same rule 
applies to p' c . 

Solving Eq. (|BT)) for /(p; P) subject to the boundary 
condition that /(p, 0) = 1 we find: 



/ 2 (p;P) 



M /S P -e 1 (p;P)-e 2 (p;P) 



Ev 



M-2e(p c ;0) 



(68) 



This completes the proof that the commutator relation 
[K, H] = iP is satisfied exactly in the subspace of eigen- 
value Ep when spin effects are omitted. 

When spin effects are included, we find that if the cm. 
frame interaction is expressed in terms of a matrix ele- 
ment involving Dirac spinors depending upon p c and p c , 
then the commutator term involving Kp„ and the spin 
terms generated by B op acting on the Dirac spinors in 
Eq. J5l| cancel. That is, let 

Vc(p'cPc) = Ul{p' c )u 2 {p' c ) V{p' c ,p c ) Ui(p c )u 2 (Pc) (69) 

where v(p' c ,p c ) is a rotationally invariant function of its 
arguments that may involve local Dirac matrix struc- 
tures, for example the Fermi covariants. The Dirac 
spinors u\, u 2l etc. in (|69|l omit the usual Pauli spinor 
factors, i.e., they are boost operators for Dirac spinors 
similar in form to Eq. (|83[1 below. It follows that the 
spinor matrix element in (|69|l involves products of terms 
like <7i • p c and a 2 ■ p c from the spinors and v. Now con- 
sider the Dirac matrix structure FiT 2 that arises from 
these products. The cancellation we want will occur, 
provided that 



B„ p r 1 r 2 -i[Kp <7 ,rir2] = o. 



(70) 



But the Dirac structures FiT 2 either involve u\ ■ a 2 or — 
as in the case of irNN coupling — factors of a n ■ p c . A 
straightforward calculation then shows that, 



and that a similar identity holds also for argument p' c 
and spin operator a 2 . Since Kp,,, commutes with u\ • er 2 , 
it follows that: 

RopV c (p'c, Pc) - i[Kp, ff , v c (p' c , pc)} = 0. 



(72) 

The remainder of the boost operator can be under- 
stood as requiring a rotation of spins when P^O, 

iP ■ (ji x p , , 

^ Wl + 2^T^ + -"' (73) 

where i = 1 or 2 and omitted terms are of order 
(66) 1/m 4 . 42j That is, the interaction in a frame with 
nonzero total momentum should have the final form 



(p';PMp;P) = /(p'; P) J R' 1J R^ 1 (p c )S 2 (p c ) 

Hp'c p c )u 1 (p c ) U2 (p c )i? 1 i? 2 /(p; P). (74) 

Factors and Ri allow 1 / m 2 parts of the terms involv- 
ing h' a and k CT to be canceled by new terms generated 
when B op acts on the spin rotation factors in Eq. Q54p. 
However, there are order 1/m errors that cannot be can- 
celed: we accept these as error terms of the boost related 
to the fact that our boost generator is not exact. 

In summary, we find that if an arbitrary cm. frame 
interaction that is rotationally and translationally invari- 
ant is constructed in other frames according to Eq. I|74|l . 
then a bound state of mass M defined by solution of the 
cm. frame equation, 

[M-2 C (p';0)]* c (p / ) = J 7&«c(p',p)*c(p), (75) 

corresponds in another frame to a state of energy Ep, 
given by solving the eigenvalue equation: 

[Bp-e 1 (p';P)-e 2 (p / ;P)]*(p / ;P) = 
d 3 p 



(2tt) e 



(p',PMp;P)*(p;P). (76) 



Wave functions for the initial and final state are now 
obtained from Eq. (|76fl using the interaction defined as in 
Eq. H74|) ■ This provides wave functions of energy Ep that 
correspond to a mass M in the cm. frame of the two- 
nucleons, and which can be inserted into Eq. (|36fl . The 
spin rotation factors are given approximately by Eq. 17311 , 
within errors of order 1/m 4 . We have omitted spin ro- 
tation effects in our calculations but we find the correct 
relation between energy and mass with high accuracy. 
The boost rule just stated is approximate because we 
have omitted terms that would be required in order to 
achieve an exact solution of the Poincare algebra. It is 
designed to be consistent with the currents derived using 
the ET reduction. 



VI. PSEUDOVECTOR-PION CONTRIBUTION 
TO THE ISOSCALAR TWO-BODY FOUR 
CURRENT 



B op (Ti • p c - i[Kp >a ,ai ■ p c ] = 0, 



(71) 



The derivation presented in Sec. II VI involves a straight- 
forward reduction of conserved four-dimensional matrix 
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elements of the electromagnetic current. It neglects re- 
tardation effects (known to be small) but also neglects the 
dependence on the relative energy po of meson-nucleon 
vertex functions. 

As we will show here, keeping the po-dependence of 
the itNN vertex results in an additional contribution 
to the 3D current that can be interpreted as a "meson- 
exchange current" contribution to the current operator. 
As stressed by Friar [3(| and Adam and Arenhovel 0] , it 
is critical that these MEC contributions be evaluated in a 
manner consistent with that used to evaluate the NN in- 
teraction. Here we implement this consistency by recall- 
ing that the instant NN potential K-et is derived by com- 
puting the ET matrix elements using a four-dimensional 
kernel and then taking the static limit \E — e — e'\ <C u> 
in the denominator. Thus in order to obtain a consistent 
Afj, for pseudovector pion coupling we should implement 
Eq. with the 4D K given by Eq. and also take 
the static limit there. This leads to a pion-range contri- 
bution to the isoscalar four current, one first identified 
by Riska a number of years ago. This contribution was 
not included in calculations of Ref. 0] and its omission 
is part of the reason that the data for A is underpre- 
dicted there. It is depicted on the right-hand side of the 
upper line of Fig. [21 For pseudovector irN coupling, the 
extra current is derived from Eq. 1(28(1 by evaluating the 
expressions using the positive-energy part of each Dirac 
propagator and the 7r-exchange kernel of Eq. (|22() . Fol- 
lowing the procedure outlined above Eq. I|24|l leads to an 
expression very similar to Eq. (|29|l , as follows 



A* — ^'et{(G'oJ^Go)K^ (G ) - 
— (G' J Pi Go)}Tet- 



(77) 



where 



^ET ~ "ET — 

9*n ■ r 2 



A^2 { 



7? 



rA - £2 



c 2 50 



2m 



"72 72 (UlU 2 , 



(78) 



with ta = \{E-p — 2e 2 + e 2 — e^). The factors in curly 
brackets in Eq. 1(78(1 are similar to those in Eq. 1)24(1 . The 
replacement of p\ — p' x by —ta occurs as the result of 
performing the integrations over time-components of mo- 
menta as dictated by the reduction to three dimensions. 
Conversely p\ — p' 2 ° gets replaced by +ta- 

The extra terms that arise due to the retention of 
p' Q — pq pieces in the numerator of the pion exchange can 
now be determined by comparing Eq. (|77|l with Eq. 1)29(1. 
This leads to the following expression, which is what we 
calculate below in order to include the PV-contact MEC: 

+(G' Q J»G Q ) (k& - K£ T ) (Go) 
+(G' Q ) (K%£-K%r) (G Q J^G )}l 



ET- 



(79) 





FIG. 2: On the left of the upper line we depict the positive- 
energy-state impulse-approximation mechanism for electron- 
deuteron scattering. On the right of that line is the PV- 
contact MEC. In the bottom line we show the p-n'y MEC, 
with the p-meson represented by the wavy line. In each case 
the blobs are the deuteron vertex functions Y. 



In fact, the dominant contribution here is due to Jo- 
Contributions to J are suppressed by p/M and are nu- 
merically small, but they are included automatically in 
our formalism. The resultant expression for the two-body 
piece of Jo agrees with Adam et al. [3l| at leading order 
in p/ra. We have also checked that the expression 1(79(1 
for this PV-coupling current agrees numerically with the 
more conventional expression as a two-body operator — 
again, to leading order in pjra. 



VII. THE p7T 7 EXCHANGE CURRENT 

Because of the quantum numbers of the deuteron, the 
p7T7 meson-exchange current is generally thought to be 
the lowest-mass mesonic excitation that makes a con- 
tribution to the electromagnetic deuteron current. The 
/97T7 exchange current can be used to provide a resonance- 
saturation model of the two-body contribution to J. This 
meson-exchange piece of Gm is of 0(eV A ) , and has an un- 

(3) 

known coupling from C^j^. The value of this unknown 
coupling can be modeled by assuming it is dominated by 
the p excitation, see e.g. Ref. 0- 

The Lagrangian governing the pnj vertex is: 



2m 



while the pNN vertex has the form 



C 



pNN 



T + i 



g p 2m 



N. 



(80) 



(81) 



This yields the two-body current depicted on the lower 
line of Fig. |2 Note that this current is automatically con- 
served. In calculating it we use the same vertex functions 
and form factors employed at the TtN and pN vertices in 
calculating the NN potential. The coupling g p7r ^ — 0.563 
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is extracted from the decay p — > This value was 

used by Hummel and Tjon |l9||. however, Truhlik, Sme- 
jkal and Khanna [^3 recently have calculated somewhat 
higher values, namely .585 ( based on p^ decay) and .610 
( based on p° decay). Although the overall sign of g p7ri 
is not determined by the experimental data, a recent lat- 
tice QCD calculation gives a posi tive value |4^. Quark 
models such as the one of Ref. |47J also provide a positive 
value for <7 p7r7 . We have used g p ^ 1 =.563 in our calcu- 
lations. For consistency with the Bonn-B potential the 
ratio fp/ g p was initially chosen to be 6.1. However, f p / g p 
is not well determined by fitting the NN data and recent 
work has suggested that a smaller value f p /g p ~ 4.5 is 
acceptable |48j . This would be closer to the f P /g P ~ 3.7 
value obtained when vector-meson dominance is used to 
explain the isovector anomalous magnetic moment of the 
nucleon We note that in nonrelativistic calcula- 

tions the f p / g p term often is neglected altogether on the 
grounds that it is higher order in 1/m |jg. In contrast, 
Hummel and Tjon performed relativistic calculations of 
deuteron form factors using f p /g p — 6.8. They found 
that the contribution of the tensor pN coupling to Gm 
cancels with the contribution of the vector pN coupling 
when |q| is about 1 GeV. This throws into question the 
validity of neglecting the tensor piece of the pNN cou- 
pling if f p /g p is as large as this. 

Even were f P /g P known precisely there is still an- 
other significant source of uncertainty in evaluating this 
exchange current: its contribution to electron-deuteron 
scattering depends crucially on the behavior of the cur- 
rent operator as a function of Q 2 , and hence on the pir^f 
form factor. In the work of Hummel and Tjon vector- 
meson dominance was used to obtain a pix^j form factor 
given solely by the u) meson: 



1 



(82) 



This same pir^f form factor is also employed in the non- 
relativistic calculations of Refs. [alii. Other calculations 
have used form factors based on quark models [2(3, • 
Such form factors tend to reduce the contribution of this 
MEC, which is also very sensitive to the cutoff masses in 
the irN and pN vertices. 

No experimental data are available to constrain the 
form factor that should be used for the pir~/ vertex. Lat- 
tice QCD calculations can determine the form factor for 
on-shell tt and p mesons. Recent lattice calculations of 
Edwards 0] indicate that for Q 2 up to 0.6 GeV 2 the 
pit"/ for factor agrees with the vector-meson dominance 
form (JS2J}. 



VIII. Z-GRAPHS 

In the full ET analysis the interactions and propaga- 
tors are defined on a complete set of Dirac plane-wave 
states: ++, H — , — h and . The bound-state equa- 
tion is solved with all possible couplings between these 



states, and the current matrix elements allow for all pos- 
sible transitions between the states. Physically, coupling 
to negative-energy states allows positive-energy, plane- 
wave solutions of the free Dirac equation to be distorted 
in the presence of a meson field. The simplest example of 
this effect — and one of the motivations for including it — 
is the shift of the mass of a Dirac particle in a uniform 
external scalar field: to — > m* = to + S. Positive-energy 
solutions of the Dirac equation in this scalar field have 
the form 



u(p) 



1 



(83) 



To recapture the effect of the replacement to — ► m* order- 
by-order in S the coupling of positive-energy states of 
mass to to negative-energy states must be included in 
the calculation. So allowing for couplings between all 
positive- and negative-energy components — as we do in 
the full ET formalism — is a general way to incorporate 
this kind of distortion of the Dirac spinors in the theory. 

It is interesting to explore the extent to which these ef- 
fects of coupling to negative-energy states are captured if 
we include them only via first-order perturbation theory 
in the nucleon- nucleon potential. Thus, instead of solving 

(|15fl in all p-spin sectors, pi,p2 = ++- H — , , we 

solve it in ++ states alone to obtain and then gen- 
erate the couplings to negative-energy states by treating 
the difference between the full interaction and the ++ 
states-only interaction in first-order perturbation theory. 
This leads to the following expression for the current ma- 
trix element: 



-Az )( u - r E T + {(G'oJ A 'G'o) 

+(G' J P G } [Ket 
+ (G'o) (K ET 



K, 



+ + ,+H 
ET 



)(Go) 
(G^G )}r++.(84) 



The second and third terms in Eq. I|84|) provide the 
leading-order corrections due to negative-energy compo- 
nents in the initial and final states, respectively. The 
leading-order Z-graph calculations we present below are 
performed in this way. 



IX. RESULTS 

With all the theoretical pieces of the puzzle assem- 
bled we now calculate electron-deuteron scattering ob- 
servables. The vertex functions employed are the ones 
calculated with all positive and negative-energy states in- 
cluded, as described in Ref. 18]. If the negative-energy 
states are dropped the interaction is exactly the Bonn-B 
potential for the Thompson equation, as derived and fit- 
ted to NN phase shifts in Ref. J25|. This model gives a 
reasonable fit to the NN data, and good deuteron static 
properties, although it is not as good a fit as some more 
recent NN potentials H,I13,0|- When negative-energy 
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states are included the deuteron binding energy changes 
slightly. To compensate for this we adjust the cr-meson 

coupling from the value of the fit in Ref. J25|, |^ = 8.08, 
to g = 8.55. 

The single-nucleon form factors are taken from the re- 
cent work of Kelly that includes data from Jefferson 
Laboratory on the ratio Ge p /Gm p and Gb„. However, 
we find that the results are very close to those obtained 
using the nucleon form factors of Mergell et al. [53] . Both 
parameterizations incorporate constraints on the asymp- 
totic shape of F\ and F% using arguments from pertur- 
bative QCD. 

We see clearly in Fig. [21 that the baseline ETIA cal- 
culation underpredicts the A data for Q 2 = 25-75 fm -2 , 
i.e. Q 2 —1-3 GeV 2 . This is true for both the curve 
labeled ET(++) — which includes only ++ states — and 
the one labeled ET(neg) — which includes the effects of 
all positive- and negative-energy states to all orders in 
Ket- Data from the two JLab experiments that have 
measured A 0, Q are denoted by triangles and squares. 
Note that these experiments confirm the trend of the 
SLAC data of Arnold et al. The extant experimental 
data for B [H^, [5(| [§3 (which, as yet, include no JLab 
data) are even less well-described. Already at Q 2 < 25 
fm -2 w 1 GeV 2 there is significant disagreement between 
our ETIA calculation and the data. 

In Fig. [31 we also present results for observables A 
and B showing the effects of the approximate dynam- 
ical boost. The curves are labeled Boost (++) and 
Boost (neg). The net effect of the dynamical boost is sig- 
nificant but not large. For either the ++ calculation or 
the one including all states, the boost shifts the minimum 
in the calculated B observable to somewhat higher Q 2 . 
As seen in Fig. [SJ the A observable is not much altered 
by the boost. These results suggest that further refine- 
ment of the approximate boost developed here, while the- 
oretically interesting, would have little phenomenological 
impact on these observables. 

A similar result is seen for observable £20 m Fig- SI 
Although effects on observables are not large, inclusion 
of the approximate dynamical boost in our calculations 
is an important step forward because it ensures that the 
deuteron mass is independent of frame, whereas an ad 
hoc factor A(P 2 ) in Ket was needed for this purpose in 
prior work. 

Calculations with and without the PV-coupling MEC 
are shown in Fig. [SI for A and B and in Fig. [SJfor £20- 
This MEC, which must be included because we chose to 
employ pseudovector irN coupling, provides a significant 
increase in the A observable, small changes in B and a 
significant shift of £20 toward smaller Q. The increase in 
A helps to move the result into better agreement with the 
experimental data. But our results for B still decrease 
too fast with Q and the minimum occurs at too low a 
value of Q. The shift of t 2 o aligns our results quite well 
with the JLab data of Abbott et al. 0,0- 

In Fig. [7| we present our results for A and B when the 




Q 2 (fm 2 ) 



FIG. 3: The deuteron structure functions A and B. The ex- 
perimental data for A of Refs. |5^.j 5 ^| are denoted by the 
open circles, while those of Refs. Q, [J are represented 
by triangles and squares. The experimental data for B of 
Refs. 55, 56] arc shown by open circles, while those of Ref. |5^ | 
are shown by squares. The dotted and short-dashed lines 
labeled ET(++) and ET(neg) are based on the ET calcu- 
lation without the dynamical boost. The long-dashed and 
solid curves labeled Boost(++) and Boost(neg) include the 
dynamical boost. 



1.0 





100.0 



Q 2 (fm -2 ) 



FIG. 4: The tensor-polarization observable £20- The older 
experimental data [58(1 and the NIKHEF data of Bouwhuis et 
al. |5fl| are shown by circles, the JLab data of Abbott et al. [Dl 
by squares and the Novosibirsk data of Nikolenko et al. |6d| 
by triangles. The lines have the same meaning as in Fig. [21 



/97T7 MEC is included using the vector-meson-dominance 
form factor (|82|l as in the work of Hummel and Tjon. 
Results are shown for three values of the tensor coupling 
ratio f p /g p = 0, 3, and 6.1 and all calculations include 
the full set of positive- and negative-energy states and the 
PV-coupling MEC. At larger values of Q, the calculated 
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0.0 20.0 40.0 60.0 80.0 100.0 

Q 2 (fm~ 2 ) 



FIG. 5: The deuteron structure functions A and B with and 
without the PV-coupling MEC. Experimental data points are 
as described in caption to Figure El 



1.0 




Q 2 (fm~ 2 ) 



FIG. 6: The tensor-polarization observable t2o with and 
without the PV-coupling MEC. Both calculations include 
all positive- and negative-energy sectors. Experimental data 
points are as described in caption of Figure 0] 



values of A are increased by the pir^f MEC. For the ratio 
f /g = 0, the /?7T7 contribution provides good agreement 
with the data for observable A. 

Meanwhile, Fig. |H1 shows results for the tensor observ- 
able £20 j which is rather insensitive to the pit"/ contri- 
bution. In contrast, this exchange current has a signif- 
icant impact on the observable B. In particular, differ- 
ent choices for the tensor pN coupling change the results 
markedly, as is evident in the lower curves of Fig. [7| It 
is clear from these results that the minimum of B is es- 
pecially sensitive to the value of f P /g P , and that experi- 
mental data are described best with zero tensor coupling. 
However, a softer pir^f form factor than that of Eq. i|82|) 




Q 2 (frrf 2 ) 



FIG. 7: The deuteron structure functions A and B. Each line 
is based on including all positive and negative-energy states, 
the boost, the PV-coupling MEC and the piry MEC. Lines 
differ only by the value of tensor pN coupling f P /g P used in 
the pivy graph. 



1.0 




0.0 25.0 50.0 75.0 100.0 

Q 2 (frrf 2 ) 



FIG. 8: The tensor-polarization observable £20- Each line is 
based on including all positive and negative-energy states, the 
boost, the PV-coupling MEC and the pn-y MEC. Lines differ 
only by the value of tensor pN coupling f /g used in the pTvy 
graph. 



would reduce the impact of the piri MEC on all observ- 
ables. Also, Hummel and Tjon [19j have estimated the ef- 
fects of a possible cuaj meson-exchange current and found 
that it can have a significant effect on the minimum of 
the B observable when a vector-meson-dominance u>a"f 
form factor is used. Clearly, a better understanding of 
the effects of meson-exchange corrections on the electron- 
deuteron B structure function is needed. Lattice QCD 
calculations along the lines of Ref. ^(| are needed for the 
pir-f form factor with Q 2 values up to 2 GeV 2 . 
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Negative-energy states have significant effects on A and 
(particularly) B. See, e.g. the curves labeled Boost(++) 
and Boost(neg) in Fig. We do not, however, find as 
large an effect from negative-energy states as was ob- 
tained in the work of van Orden, Devine and Gross [2(J . 
There the inclusion of such effects produced an accept- 
able description of the B data. This is likely due in 
part to the fact that we use a pseudovector irN coupling, 
whereas in Ref. [2(j an admixture of about 25% pseu- 
doscalar irN coupling was employed. Pseudoscalar irN 
coupling is known to produce larger effects from negative- 
energy states. 

In order to clarify the significance of including all pos- 
itive and negative-energy states in our analysis, we also 
calculate the leading-order Z-graphs. Figs. El and 117)1 
compare results obtained using ++ states only, adding 
the leading-order Z-graphs as in Eq. Q84[). and including 
all positive and negative-energy states. Meson-exchange 
currents are omitted. In the case of observable B the 
leading-order Z-graphs do not accurately reproduce the 
results of a calculation which includes all positive- and 
negative-energy states to all orders in perturbation the- 
ory. First-order perturbation theory for Z-graphs is ade- 
quate for observables A and £20, but not for B. Solving 
for the deuteron vertex functions non-perturbatively in 
all positive- and negative-energy sectors one finds smaller 
effects on the B observable due to couplings to negative- 
energy states than are predicted by leading-order Z- 
graphs. 



the PV-coupling current are implemented in a consis- 
tent formalism and calculation. However, the data for 
the B observable are not explained unless there is some 
non-standard contribution, such as the pirj MEC with 
f p = 0. Other analyses have suggested that various rela- 
tivistic effects, effects of negative-energy states or effects 
of the PV-coupling current might be larger than we find, 
but we believe our examination of each of these effects to 
be reliable within the ET formalism. 



If a vector-meson-dominance pirj form factor is cor- 
rect, then the pirj MEC can help to explain the B data 
provided that a very small tensor pN coupling is used. 
The ratio f p /g p should be substantially less than the 6.1 
used in one-boson-exchange models of the NN interac- 
tion, or even a value which would explain the isovector 
magnetic moment of the nucleon. Fig. [S] shows that the 
t2o data of Refs. Q, |^3, EH |(5(j are well-described by our 
approach out to Q 2 ~ 2GeV 2 . This observable is fairly 
insensitive to some of the dynamics that plays a role in 
A and B (single-nucleon form factors, the pnj MEC). 
However, it is quite sensitive to relativistic effects (see, 
for instance H,E3): so it is gratifying that our approach 
reproduces the data, especially that of Ref. pj at large 
Q 2 , so well. 
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FIG. 9: The deuteron structure functions A and B. In all 
cases the boost is included and the meson-exchange correc- 
tions are omitted. Dotted lines are based on including ++ 
states only, solid lines are based on including all positive and 
negative-energy states and dashed lines are based on including 
++ states and adding the Z graphs in first-order perturbation 
theory as in Eq. 1841 . 

Our results suggest that the existing data for the 
deuteron's A and £20 observables are described reason- 
ably well when relativistic dynamics, boost effects and 
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